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Abstract 

We give an expression for the partition function of a one-dimensional log-gas comprised of 
particles of (possibly) different integer charge at inverse temperature [} = 1 (restricted to the line 
in the presence of a neutralizing field) in terms of the Berezin integral of an associated non- 
homogeneous alternating tensor. This is the analog of the de Bruijn integral identities [3] (for 
y6 = 1 and ft = 4) ensembles extended to multicomponent ensembles. 
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1 Introduction 

We imagine a finite number of charged particles interacting logarithmically on an infinite wire mod- 
elled by the real line. Different particles may have different charges (which we will assume are 
positive integers), but any two particles with the same charge, that is, of the same species, are as- 
sumed to be indistinguishable. A potential is placed on the wire to keep the particles from escaping 
to infinity. This system is placed in contact with a heat reservoir with inverse temperature (3. 
We will consider two ensembles: 

1 . The Canonical Ensemble. We assume that the number of each species of particle is fixed. 

2. The Grand Canonical Ensemble. We assume that the sum of the charges, that is the total 
charge of the system, is fixed but the number of each species is variable. 1 

Our goal is to provide a closed form of the partition function of these ensembles, for certain values 
of /3, in terms of Berezin integrals. As is standard, we will find that the partition function for the 
Grand Canonical Ensemble is the generating function for the Canonical ensemble as a function of 
fugacities of the species of particles. After a minor modification, the partition function can also be 
seen as the generating function for the correlation functions of both the Canonical Ensemble and the 
Grand Canonical Ensemble. 



'This research was supported in part by the National Science Foundation (DMS-0801243) 

'The standard notion of the Grand Canonical Ensemble is that where the number of particles is not fixed. That is, in its 
traditional sense, the Grand Canonical Ensemble is the direct sum over all possible values of the sum of the charges. What 
we refer to as the Grand Canonical Ensemble might be better referred to as an isocharge or zero current Grand Canonical 
Ensemble. 
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1.1 The Setup 

Let J > be an integer and suppose q = (q\,q2, ■ ■ ■ , <?./) is a vector of positive integer (charges) with 
each of the qj distinct. We imagine a system of particles consisting of M\ indistinguishable particles 
of charge q\, M 2 indistinguishable particles of charge q 2 and so on. We will refer to q as the charge 
vector and M = (Mi, M 2 , . . . , Mj) as the population vector of the system. 

These particles are restricted to lie on an infinite wire, identified with the real axis 2 , and interact 
logarithmically, so that the energy contributed to the system by a pair particles with charges q and 
q' located at x and x' is given by-qq' log \x - x'\. (Infinite energy is allowed in the situation where 
x = x'). We suppose that the particles of charge q\ are identified with the location vector x 1 = 
(x\,x\,..., x l Mi ); the location vectors x 2 , . . . , x y are similarly defined. If Mj = for some j then x J 
is taken to be the empty vector. The particles are placed in a neutralizing field with potential U so 
that the total potential energy of the system is given by 

J Mj J _ Mj M t 

Euix 1 , x 2 , . . . , x J ) = qj 2 U{x ] m ) - Yj <r) J] 1o 8 |x " ~ x '»! ~Yj qjqk YjYj l0g |X " ~ 

j=\ m=l j=l m<n j<k m=l n=l 

We assume that the system is in contact with a heat reservoir at inverse temperature /3, but energy 
is allowed to flow between the reservoir and the system of particles. In this situation the Boltzmann 
factor, which gives the relative density of states, is given by 

Q M (x\x 2 , . . . ,x y ) = e-^x'.-.V) (L1) 

J Mj J _ Mj M k 

= { n n ^ qju(4i) } x { n n ^ - ^} x { n n n & - 

j=l m=l j=\ m<n j<k m=l n=\ 

The probability (density) of finding the system in a state determined by the location vectors x 1 , x 2 , 
. . ., x J is then given by 

! 2 j Q M (X 1 ,X 2 ,...,X / ) 



/? M (x 1 ,x z ,...,xO = 



Z M M x \M 2 \---MjV 
where the partition function of the system is given by 



Z m = 777777 777 f "T ^(x 1 ^ 2 , . . . ,x J )dn M ^)dp. M ^ 2 ) ■ ■ ■ dfi M '(x J ), (1.2) 

Mi!M 2 ! • --M/! J r m, J r mj 

and p M is Lebesgue measure on R M . 3 The factors of M\\M 2 \ •••M J \ appear since a generic state 
of the system x 1 , x 2 , . . . , x J has this many different representatives. We will always assume that the 
external potential U is such that Z M is finite. 

For the grand canonical ensemble, We may view M as a random vector and the probability 
(density) of finding the system with a prescribed population vector M and state x\x 2 , . . . ,x J is 



2 With a minor modification, much of what is presented here can be shown mutadis mutandis for multicomponent log- 
gasses confined to the unit circle. See [5] for the circular case for /} = 2, and [6] for some physical application. 



3 If Mj = for some j then we will use the convention that 

n M (x' , x 2 , . . . , x J ) dn M i(xi) = n M (x' , x 2 , . . . , x'); 



X* 



alternately, in this situation we may assume that the integral over R M ' does not actually appear in our expression for Zm- 
Likewise we will assume that sums and products over empty sets are respectively taken to be and 1 . 
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given by /?m(x',x 2 . . . ,x J ) ■ prob(M). Classically, the probability of finding the system in a state 
with (allowed) population vector M is taken to be 

P rob(M)=zf'zf ...zf 7 ^, (1.3) 

where 

Z N = 2_, z l z 2 ...Zj Z M 

M 
Mq=» 

and z = (zi , . . . , zj) is a vector of positive real numbers called the fugacity vector. 

It shall sometimes be convenient to view z as a vector of indeterminants and Z# = Z#(z) as a 
polynomial in these indeterminants. Our main result will be to show that, for certain values of ft 
and q, Z N (z) can be expressed as a Berezin integral with respect to the volume form in R N of the 
exponential of an (explicitly given) alternating element (i.e. form) in the exterior algebra A(R. N ). By 
construction, Zm is the coefficient of zf' ■■ ■ zf' , and thus the integral formulation of Z N (z) is exactly 
the generating function we seek. 



2 Wronskians, Berezin Integrals and Hyperpfaffians 

Here we collect the machinery necessary to state our main results. 

Given a non-negative integer L, let L = { 1 , 2, . . . L], and, assuming K > L is an integer, let 

be a strictly increasing function, 

< t(l) < t(2) < ••• < t(L) < K. 

We will use such functions to keep track of minors of matrices, elements in exterior algebras and 
Wronskians of families of polynomials. Such indexing functions will always be written as fraktur 
minuscules. 



2.1 Wronskians 

A complete family of monic polynomials is a sequence of polynomials 3> = (p\,p2, ■ ■ •) such that 
each p„ is monic and deg p„ = n - 1. We define the L-tuple % = (pt(i), ■ ■ - ,pt(L))- And, given 
< £ < L we define the modified fth differentiation operator by 

D°f(x)=f(x) and D e f(x) = (2.1) 
The Wronskian of T t is then defined to be 

m{%-,x) = tet[D i - 1 p m (x)\ k{=i . 

The Wronskian is often defined without the l\ in the denominator of (2.1); this combinatorial factor 
will prove convenient in the sequel. The reader has likely seen Wronskians in elementary differential 
equations, where they are used to test for linear dependence of solutions. 
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2.2 The Berezin Integral 

If d , . . . , Ck is a basis for R K , then et = et(i) A • • • A et(L) is an element in A L (R K ), and {et 1 1 : L f K) 
is a basis for A L (R^). In particular, we will denote 

e vo i = ei A e 2 A • • • A e^. 

Given < k < K we define the linear operator d/de k : A L (R K ) -> A^^R^) by 

q < (-l) ff e t( i) A • • • A e t ( a -i) A e t(a+ i) A • • • A e t( L) if k = r l (a); 

9ek \ otherwise. 

That is, if appears in $ then det/de k is formed by shuffling to the front of e t (taking into account 
the alternation of signs) and then dropping it. Given < k\, . . . , &m ^ K we then define the Berezin 
integral as the linear operator on A(R K ) — > A(R Ar ) specified by 



e t de ki de k2 ■ ■ ■ de ku = — - — e t . 

de kM de k2 de kl 



Berezin integrals were introduced in [1] as a Fermionic analog to the Gaussian integrals which 
appear in Bosonic field theory. 

We will mostly be interested in Berezin integrals of the form 



J e t de mi = J 



et dei ■ ■ ■ d^K- 



In this case, the Berezin integral is simply the projection operator K{R K ) — > A K (R. K ) = R. Notice in 
particular that, if <x e Sk then 



e<r(i) A A • • • A e^K) de m \ = sgn <r. 



2.3 Exponentials of Forms and Hyperpfafflans 

Given cj e A(R K ) we define w A ° = 1 and for < m 



Using this we define 

°° . .Am 

e" = y 

m=0 

If a) = ojq + a>\ H + ojk with a> k e A k (R K ) then it is easily verified that 

e " = e"° Ae U| A • • • A e" K . 



Moreover, e"° is a real number equal to its traditional definition, and if k > then the sum defining 
e" k is a finite sum. 

In the situation where k divides K, that is K = km, then we define the hyperpfaffian PF(aj k ) to be 
the real number defined by 

= PF(w*) e vo i. 

ml 
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Alternately, 

PF(^) = j ef* e vol . 

The hyperpfaffian is related to the Pfaffian of an antisymmetric K x K matrix by associating 
the matrix to a 2-form in the obvious manner. We see therefore that the Berezin integral formed 
with respect to e vo i is a generalization of hyperpfaffians, which themselves are generalizations of 
Pfaffians. 

3 Statement of Results 

Suppose b is a positive integer and f3 = b 2 . Set K = bN and Lj = bqj for j = 1,2,...,J. For any 
complete family of monic polynomials T we define oj\, on, ■ ■ ■ , 0)j e A(R K ) as follows. 

1. If Lj is even, 

coj = Y ( f e -^WWr(T t ;x)rfxL; (3.1) 

2. If Lj is odd, 

ojj = V (i f r e -^ J/w e-^ J/ WWr(Tt;x)Wr(J' u ;y)sgn(y-x)c/xd/y}et A e u . (3.2) 

Notice that coj is in A L '(R K ) when Lj is even and is in A 2L i(R. K ) when Lj is odd. 
Theorem 3.1. Suppose f5 — b 2 and K = bN is even. Given a charge vector q let 

Lj = bqj; 7 = 1,2,...,/, 
and, for any complete family of monic polynomials, define the form cj e A(R hN ) by 



Xz) = J^ZjOJj, 



7=1 

where ojj is defined as in (3.1) or (3.2). If the Lj are positive integers, at most one of which is odd, 
then 

Z N (z)= f e^d€ ml . 



Remark. This is an algebraic identity which can be written more generally by replacing the integral 
over R with integrals over other sets (for instance, the partition functions for multicharge circular 
ensembles can be likewise expressed in terms of Berezin integrals). The only analytic prerequisite 
is the finiteness of the Zm which allows for the use of Fubini's Theorem. 

This theorem covers certain situations which have appeared before. Certainly the Pfaffian parti- 
tion functions of the classical one-species ensembles GOE and GSE (and their non-Gaussian vari- 
ants) are a corollary. These cases follow from the 'classical' de Bruijn identities [3]. (see [7] and the 
references contained therein for their applications to random matrix theory). The classical ensem- 
bles can be viewed either as ensembles of charge 1 particles at respective inverse temperatures (3 = 1 
and f3 = 4, or to ensembles of charge 1 and charge 2 particles (respectively) at inverse temperature 
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P — 1. Recent work by the author, B. Rider and Y. Xu produced (among other things) a Pfaffian 
formulation of the partition function for the grand canonical ensemble for charge vector q = (1,2) 
at inverse temperature ft = 1 [9]. When the external field is Gaussian, and by tuning the fugacity, 
this latter ensemble provides an unusual interpolation between the classical ensembles GOE and 
GSE. Moreover, the skew-orthogonal polynomials necessary to solve the ensemble (that is explic- 
itly derive the matrix kernel in which the correlation functions can be expressed and analyzed) were 
explicitly computed in terms of certain generalized Laguerre polynomials. Amongst other results, 
this allowed us to compute the distribution of the number of each type of particle for various fugac- 
ities. This analysis follows similar work for the two charge circular ensemble with charge vector 
q = (1, 2) initiated by R Forrester and others (see [4, §7.10] and the references therein), and indeed 
the partition functions for those ensembles can be expressed as a Pfaffian, and hence in terms of 
(variants) of the above Berezin integrals. 

Recent work of the author [10] has lead to a hyperpfaffian expression for the partition functions 
of single-species ensembles of charge 1 particles when [3 = L 2 is a perfect square, or ft = L 2 + 1 is 
even. In the former case, these ensembles can also be interpreted as systems of charge L particles at 
/J=l. 



3.1 Correlation Functions 

Using a slight modification, the partition function gives a generating function for the correlation 
functions. For single species ensembles, the correlation functions are simply renormalized marginal 
densities. For multicomponent ensembles, however, the situation is more complicated (though the 
marginal probabilities are an important ingredient). 

For fixed population vector M and vector m = {m\,m2, . ■ . ,mj) with < mj < M, we define 

r = (£',...,£,) and y J = (y J v .-.,y j Mnm ), 

and set 

=({{,..., {L r yi,...,y j Mi _ mj ). 
The mth marginal probability density of pm is then given by 

PmA%\t,--;% J )= f "■ f PM(l l Vy\...,l J Vy J )d t i M >- m >(y 1 )---d t i M >- m i(y J ), 

and by symmetry, the probability (density) that our system is in a state (x 1 , . . . ,x J ) which occupies 
the substate (|j, . . . , | y ) (that is, viewed as sets, | ; c x j for each j) is given by 

^M,m(§ , | ,...,?) = — • • • — 77/>M,m(? ,£=,...,§) 

(Mi - mi)! (Mj - mj)\ 
1 



Z M (M X - mi)! • • • {Mj - mj)\ 



f ... f Qm^ 1 Vy 1 ,...,i / Vy y )^ M '-" 1 '(y 1 )---^ m/ (y y )- (3.3) 



X 



This is the mth correlation function for the canonical ensemble with population vector M. 

To get the mth correlation function for the grand canonical ensemble we need to sum over the 
related correlation function for the canonical ensemble over all allowable population vectors M with 
nij < Mj for each j (a situation we will abbreviate by m < M), taking into account the probability 
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of being in a state with prescribed population vector. That is, the probability (density) of the (grand 
canonical) system is in a state (x 1 , . . . ,x y ) which occupies the substate (% lt ..., | y ) is given by 

probftli, . . . , £,) c (x 1 , . . . ,x y )) = J] prob(M) • R M , m &\¥, • • • - §')• 

M.q=A' 
M>m 

Denoting this density by Rn,™, (1-3) and (3.3) yield 

j z Ml • • • z Mj 

Z N( Z ) ( M l - • • • (My - OT 7 )! 

M>m 

x f ■■• f Q M (? 1 Vy 1 ,...,^Vy y )J J u Ml -'" 1 (y 1 )---^ M/ " m V)- 

Notice that, by ignoring the prefactor Z N (z), and up to an easily recoverable constant, the coefficient 
of z^ 1 • • • zf 1 in is the mth correlation function for the corresponding canonical ensemble. 

We can in turn give a generating function for the correlation functions for the grand canonical 
ensemble as follows: Let c-' = (c } v c y , . . . , c J N ) and tj = (£ y , £ y , • • • , £ 3 N ) and define the measures Vj 
and 77^ by 

d N 

= e-^M and ^.(x) = e^^W V c ^(x - £('), 

where £(x) is the probability measure with unit mass at x — 0. 

It is convenient at this point to index the forms ojj from (3.1) and (3.2) by vj so that, for instance 
when Lj is even, 



Quantities which are dependent on v = (vi, . . ., vj) will be denoted by, for instance, Zj^,Z^ and 
w v (z). Theorem 3.1 is purely algebraic, and thus, we have, for instance that 

Zl(z) = I exp{w v (z)} de mh 



We can generalize these quantities by replacing v with other vectors of measures. The following 
theorem gives particular relevance to Z^ +1| (z, c 1 , . . . , c y ), where the notation indicates the additional 
dependence on the c 7 . 

Claim 3.2. The mth correlation function of the grand canonical ensemble is the coefficient of 



W lU ZM) ■ 

7=1 e=i ZjnW 



That is, if% J = (f y , . . . , £j, ), and we define 



ffnj 



dcJ dc\ dc J mj 

then, 



Z N (z) 



dc 1 dc J 



z; + Vc',..., c y ) 



c'=- -=c'=0 
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Moreover, the mth correlation function of the canonical ensemble with population vector M is given 
by 



fllUmG 1 ,...,^) 



1 



Z M M X \ ■ ■ ■ Mjl dzf 1 dzf'ldc 1 dc' 



gmi gmj 



zrvc 1 ,...,^) 



z=c'= -=c'=0 



The proof of this claim is standard (it is the multicomponent version of the 'functional differen- 
tiation' method), and follows mutatis mutandis that for Ginibre's real ensemble [2, Prop. 6]. 

To write the correlation functions explicitly in terms of a Berezin integral (taking all of the Lj to 
be even for convenience), we note that 

o?r i= X { f Wr(^tMv; + 77>)W=^+^', 



and 



Hence, exp{w v+11 } = exp{<j/} A expfa/ 1 }. 

This is useful, since the first term in the right hand side is independent of the c- 7 . The following 
maneuvers are elementary 

ffni ffn, 

— exp{w v+ i} = exp{w v ) A — expi^/ 1 } 



5c 1 



mi 

= exp{w v ) A expW) A j f\ e-Pi< u ^ ^ Wr(3> t ; fi)e t \; 
note that since all forms are even, we do not have to specify their order. It follows that 



ffn 



— exp{w v ) A expfw* 1 : 
5c 1 



c'=0 



exp{c<j v } A 



e -PqiU(( l e ) 
1=1 V.L. 'K 



J] Wr(? t ;^)e t }, 



and that 



ffni ff"j 

_..._exp{ w v }Aexp{^; 



dc' 



c'=- "=c'=0 



= exp{w v ) A 



J nij 

e -tojui£l) 2 Wr(? t ;^)et - 

7=1 (=1 tLjSK 



We therefore have the following corollary to Claim 3.2. 

Corollary 3.3. If for each 1 < j < J, Lj is even, and % j = (£[, . . . , e ^ m '< then 



R N , m &, ■■■,¥) = ^tt f exp{o/}A 



7 m j 



,-Pqim't) 



7=1 f=l 



2 WrCP t ;$et ^voi 



Note that we do not have to justify the exchange of the derivatives and the 'integral' in Claim 3.2, 
since the Berezin integral is not an integral in the traditional sense. That is, Claim 3.2 is an algebraic, 
not an analytic, identity. Notice also, that the quantity in braces is an (m • L)-form, and therefore 
only the projection of exp of onto the space of (K - m • L)-forms will make a contribution to the 
mth correlation function. Finally, we note that a similar formula for the partial correlation function 
^M,m is available via functional differentiation with respect to the z variables. 
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4 The Proof of Theorem 3. 1 

4.1 The Confluent Vandermonde Determinant 

Suppose 0<L<^,ieR and T = {pi,pi, ■ ■ •) is any complete family of monic polynomials. We 
define the K x L matrix 

V*(z) = [oWxC, . 

and given an admissible population vector M and x 1 , . . . , x y with x J e R Mj , we define the K x K 
confluent Vandermonde matrix by 



\ M (x\...,x J ) 



V L '(x|) ••• V Ll (x l Mi ) ■■■ \ Lj (xi) ■■■ Y Lj (x j Mj ). 



Mi Mj 

(Recall that Lj = V^<?j)- In this case, the confluent Vandermonde determinant identity [8] has that 

J , , M i M k 



det V M (x', . . . = { f] Y]U - 4^} x { f[ f[ FN - x - )L ' Lk ) (4 ' 1} 

j=\ m<n j<k m=l n=\ 

= { n - ^ x { n ft n<* - ^ fmt \ 



7=1 m<n j<k m=l n=\ 

When all of the Lj are even, it follows from (1.1) that 

. / ^ 

Qm(x',x 2 , ...x'»-{| II I e-^^Jdet V M (x', . . . ,x J ) (4.2) 

j=l m=l 

We will deal with the situation where one of the Lj is odd in Section 4.5. 

4.2 The Laplace Expansion of the Determinant 

Each i ■ L y K specifies a unique t' : K - L S N_ whose range is disjoint from t. Given a K X K 
matrix V = [v m> „] and t, u : L /* K then we may create aLxL minor of V by selecting the rows and 
columns from the ranges of t and u. That is, we write 

V t , u = [vt(k)M{)]k,e=\ ■ 

Notice that the complementary minor to Vt iU is given by Vf >U '- 
We define sgn t by 

sgnt = J e t A e v de m \. 

More generally let 

t } m '■ Lj ~* K. where j = 1 , 2, ... J and m = 1,2,..., Mj, 

and set 

?=(t|,...,t Mi ,...,t(,._4) 

Mj 
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We will use t to select minors of V M (x 1 , . . . ,x J ) each of which depends only on a single location 
variable. We denote the set of all such t by 3 m 
We define sgn t by 

sgn t = j £jj A • ■ ■ A A • • • A e t ; A • • • A e t ^ de m \. (4.3) 

Mi Mj 

Clearly, sgn t = unless the ranges of the various t J m are mutually disjoint, and otherwise sgn t is the 
signature of the permutation defined by concatenating the ranges of the various t J m in the appropriate 
order. 

We will reserve the symbol i for the vector whose coordinate functions are given by 

\Ut) = t + (m - l)Lj + MiLi + • • • + Mj-\Lj-\ 

That is, for instance, if L = (2,3) and M = (2,2) then the ranges of tj.t^tj and x\ are given 
respectively by {1,2}, {3,4}, {5,6,7} and {8,9, 10). Clearly sgn? = 1. 

This notation is convenient to represent the Laplace expansion of the determinant (which we will 
write in the form most useful for our ultimate goal). 

/ M i 

det V = ? f[ \~\ det . (4.4) 

te3 M j=l m=l 

Applying (4.4) to V M (x' , . . . , x y ) we find 

J M, 

det V M (x> , . . . , x J ) = J] sgn t f] \] det Vj^xi), 

where the notation reflects the fact that V^ 1 . (x J m ) in independent of all location variables except x J m . 
From the definition of V M we see that det Vjf . ; (x J m ) = Wr( "P^ ; x J m ), and therefore 

J Mj 

det V M (x' , . . . , x y ) = Yj s 8 n ? n PI Wr(J> t,'„ ; x ->- ( 4 ' 5 > 

Fe3 H J=1 m=1 

4.3 Fubini's Theorem 

From (1.2), (4.2) and (4.5) we have that 

J Mj 



/.\\ — 

' Mi!M 2 ! 



^ sgn ? f f ■■■[ flfl e-^^Wri?, ; xjj 
\---Mj\ f-i Jrmj J r m 2 J rMj ii ii 

x dn Ml (x l )dn M \x 2 ) ■ ■ ■ d^ Mj {x J ). 

Fubini's Theorem implies then that 



i 1 Mj r 

ZM = -VM 2 !...M y ! l> r nn le-^Wr(? t r,x)dx. (4.6) 



M,. 

te3 M 7=1 m=l 
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Thus, 



Mi Mi Mi 



J Mj 



IV1 „_ K! IfcOM J 



4.4 Enter the Berezin Integral 

Using the definition of sgn t (4.3) we find 



Ziv(z) 



A • • • A &i A • • • A e t J A • • • A e t J de m \ \ 

M] ] Mj J 



7=1 ■'" m=l ^ R 



Exploiting the linearity of the Berezin integral, 



^ = f[ Z Z A ^ A^^W^;,)^ 

Mq=« ItJ M 



d&,oh 



where the wedge products are taken in the standard order. 
Next we may expand the sum over t e 3m as 



so that 



z<->- Z Z ■■■ Z (-)• 



^ = /[ZA^t Z A.(f^^ P )^ 



de m \. 



We observe that 



Mj 

V Azjl \ e-^ u(x) Wr('? i j;x)dx\e i i 



= ( z ; z ( r ^Vti^ir 

t:L.//f ^ R 



and hence 



«•-/[? w 
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Now, we can remove the restriction M • q = N from the sum in this expression, since the Berezin 
integral will be zero for any M not satisfying this condition. (If M does not satisfy this condition the 
form in the integrand will not be in A K (R K ) and hence its projection onto A K (R K ) s R will be 0). 
Thus, 

^/[zA^N 

J M j=l J' 

-/ ZZ-ZA^K' 

^ M,=0M 2 =0 M,=0 7=1 ■'' 
^ 7=1 M=l 

= e ZlWl A e Z2 " 2 A • • • A e z ^<ie vo i 
= J e tt(z W vol , 

as desired. 



4.5 When one of the Lj is odd 

In the case where exactly one of the Lj is odd, we will reorder the qj so that L\ is odd and La, ■ ■ ■ , Lj 
are even. In this situation, (1.1) and (4.1) imply that 

J Mj 

ilMlx'.x' *'>-{ | | I |<' /; "' : ' '}{ [1 sgn^-x^jdetV^x 1 ^ 2 ,...^ 7 ), 

7=1 m=l l<m<n<Mi 

where the additional factors of the form sgn(x^ - x} n ) exist in order to make the expression non- 
negative for all choices of x 1 , x 2 , . . . , x 7 . Defining the Mi x Mi antisymmetric matrix 

When K is even, so is Mi, and in this situation 

PfTXx 1 ^ Y] ^nixl-xl). 

\<m<n<M\ 

Thus, 

Q M (x 1 ,x 2 ,...,x' / ) = j P| ["[e-^^HPf^x 1 ) detV M (x 1 ,x 2 ,...,x y ), 

7=1 m=l 

Following the analysis of the case where all Lj even we find the analog of (4.6) in the current 
situation is 

1 i 3 M ' r 1 

ZM = M l! M 2 !-M /! Z sgnrfnn J/ WW ^^)^} 

x PfT(x) fle^'W^n,) UyW Ml (x) 
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And, 



Mi M 2 Mj 



J Mj 



Mi\M 2 \---Mj\ 



Mq=» 



r r Ml i 

Pf T(x) Ff e"^ 1 ^'Wrf?,, ; x m ) d^ M> (x). 



Using the same maneuvers as before, we can write 
{ {ZjlOjT M < -'" 



Zm = J z.A^tt^ L 

J ul^- 2 1 1 •! 

r Ml 

( J rM] Pf T(x){fj e -^ t/fe) Wr(0' t i ;x m )}^ Ml (x))e t| A e t i A • • • A 



t/e, 



vol- 



It is shown in [10, Section 4.2] that 

-M, . „ . Mi 



Ml V.ilV jRMl l U j y ' 



(zi^i) AMl 



(The left hand side of this expression is the partition function of a system of M\ particles each of 
charge q\ when ft is an odd square; showing partition functions of such systems is a hyperpfaffian 
was one of the goals of [10]). 
We therefore have that 



de m \ 



J(z W vol , 



as desired. 
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